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Abstract: We study Higgs production and decays in the context of natural SUSY, allowing for an 
extended Higgs sector to account for a 125 GeV lightest Higgs boson. Under broad assumptions, Higgs 
observables at the LHC depend on at most four free parameters with restricted numerical ranges. Two 
parameters suffice to describe MSSM particle loops. The MSSM loop contribution to the diphoton 
rate is constrained from above by direct stop and chargino searches and by electroweak precision tests. 
Naturalness, in particular in demanding that rare B decays remain consistent with experiment without 
fine-tuned cancellations, provides a lower (upper) bound to the stop contribution to the Higgs-gluon 
coupling (Higgs mass). Two parameters suffice to describe Higgs mixing, even in the presence of loop 
induced non-holomorphic Yukawa couplings. Generic classes of MSSM extensions, that address the 
fine-tuning problem, predict sizable modifications to the effective bottom Yukawa yb- Non-decoupling 
gauge extensions enhance yb, while a heavy SM singlet reduces yb- A factor of 4-6 enhancement in the 
diphoton rate at the LHC, compared to the SM prediction, can be accommodated. The ratio of the 
enhancements in the diphoton vs. the WW and ZZ channels cannot exceed 1.4. The h — > bb rate in 
associated production cannot exceed the SM rate by more than 50%. 



Contents 



1 Introduction 1 

2 Loop effects 3 

2.1 Stops 4 

2.1.1 Higgs mass and fine-tuning 4 

2.1.2 Higgs couplings to gluons and photons 5 

2.1.3 Large stop mixing vs. fine-tuning in BR(B — > X s j) 6 

2.1.4 Summary of stop effects 8 

2.2 Charginos 8 

2.3 Charged Higgs 10 

2.4 Staus and sbottoms 10 

3 Higgs mixing 11 

3.1 Two Higgs Doublet Model 11 

3.2 MSSM analysis 13 

3.3 Non-decoupling Z?-term models 15 

3.4 Non-decoupling _F-term models 16 

4 Natural SUSY predicts 18 

4.1 Interpreting current data 21 

5 Conclusions 21 
A More details of _D-term models 23 



1 Introduction 

Current LHC data have provided no evidence for beyond Standard Model (SM) physics, in particular, 
supersymmetry (SUSY). Instead, lower bounds were set on weak-scale SUSY, especially on the masses 
of the gluino and first two generation squarks [1, 2]. It is not inconceivable that, with time, we will 
learn that the weak scale is tuned [3]. 

Nevertheless, the game is not over for naturalness [4]. Light stops and sbottoms have only been 
excluded up to about 200 GeV [4-9]. Stronger limits exist but are model dependent, and it is still easy 
to find regions in the squark-neutralino-chargino parameter space where these limits cease to apply, 
with no need to invoke tuned degeneracies. On the theory side, rather than excluding naturalness, it is 
a consistent approach to take the experimental constraints as guidelines for model building of natural 
soft SUSY breaking. In particular, flavored SUSY breaking, with squarks of the first two generations 
much heavier than those of the third, gains motivation. Attempts to achieve this spectrum include [10- 
15]. It is also conceivable that SUSY is natural and all squark generations are light, but the collider 
signatures are altered so that the current searches are evaded. For example, the classic missing energy 
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signature may not apply. Models in this spirit include i?-parity violating SUSY [16, 17], stealth 
SUSY [18, 19], compressed SUSY [20-23] and supersoft SUSY [24]. 

In light of this debate between paradigms, the recent hint of a 125 GeV Higgs, reported by 
CMS [25-27], ATLAS [28-30] and possibly also by the Tevatron experiments [31], is exciting [32-58]. 
In the minimal supersymmetric SM (MSSM), a 125 GeV Higgs would imply fine-tuning at the few 
parts per thousand or so [59, 60]. If SUSY is to be natural then it cannot be minimal, namely, new 
interactions beyond the MSSM must deform the Higgs sector. 

Even if natural SUSY exists, given its evasive nature so far there is no guarantee that we shall 
see direct signatures of new physics any time soon. However, data have already begun to accumulate 
for one particular new set of measurements, namely, the production cross sections and decay rates 
of the Higgs boson. This data set is rich: a 125 GeV SM-like Higgs supports O(10) experimentally 
independent production and decay channels, including e.g. gluon fusion (GF) and vector boson fusion 
(VBF) production and decays such as h — > 77, WW, ZZ, bb, tt. This upcoming information pro- 
vides promising means to probe indirectly the existence and details of natural SUSY. Given that the 
experimental uncertainties for each individual channel are large, the degree by which one can draw 
information from this new data set depends crucially on the degree of predictive power provided by 
the theory. This sets the scope for our current analysis. 

In the MSSM, the third-generation sfermion, higgsino, gaugino, and Higgs sectors contain many 
free soft SUSY breaking parameters that can affect the Higgs couplings. The question arise what 
definite predictions can be made if naturalness (together with experimental constraints) is used as 
guide. The question gains depth given that a non-minimal Higgs sector would plausibly introduce 
additional free parameters in the scalar potential. 

In this paper we attempt to partially address this question. Our working assumptions are: 

1. We assume that an mh «125 GeV Higgs exists, and is embedded in a natural (and therefore 
non-minimal) supersymmetric model; 

2. While we do not commit to the MSSM quartic potential, we still assume that the weak scale 
Higgs sector is described by an approximately type-II two Higgs doublet model (2HDM); 

3. We neglect beyond-MSSM physics in loops, assuming that the leading loop corrections to the 
Higgs-SM couplings involve MSSM fields only. 

None of these assumptions is necessarily true, but with them we obtain a predictive framework and 
show that it can be definitively tested against Higgs data. 

The outline of the paper is as follows. In general, we study the modifications to the Higgs couplings 
to SM particles, denoted by 



with i = t,V, G, 7, b, r standing for top, massive vector gauge boson, gluon, photon, bottom and 
tau respectively 1 . Our main task is to organize existing computations into concrete predictions, 
distinguishing relevant from irrelevant contributions in our framework. 

In Sec. 2 we consider MSSM loops including effects from stops, sbottoms, staus, gauginos, hig- 
gsinos and Higgs. Defining a measure of fine-tuning and demanding it not to exceed about 1:10, we 
argue that only stops and charginos can induce coupling modifications larger than ~ 5%. The most 

1 While this may not be strictly necessitated by data [61], we make the further assumption that = Tz = ry. 
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relevant constraints on the stop contribution to tq are electroweak precision tests and direct searches. 
Demanding less than 4a tension in (Ap/p), the correction is limited to r G < 1.3. Rare B decays con- 
strain stop mixing and we incorporate BR{B — > X s ~/) into our assessment of the model fine-tuning, 
finding that r G > 0.85 is preferred for naturalness. Very light charginos can vary the vertex in 
the range 0.7 < r| < 1.1; the effect is limited by direct searches. 

In Sec. 3 we consider Higgs mixing. We show that once one assumes that the weak-scale Higgs 
sector is an approximately type-II 2HDM, with natural flavor conservation broken only by tan /3- 
enhanced, but otherwise small, loop effects, then the analysis of Higgs observables involves only four 
free parameters. Because of the ~ 64% branching fraction of a 125 GeV SM-like Higgs to bb, the 
coupling rb is of key phenomenological importance. We estimate the possible variations in r^, expected 
in concrete SUSY models that fall in our framework. Non-decoupling D-term models, obtained by 
integrating out additional gauge interactions, generically predict r& > 1. The precise result depends 
on the mass of the heavier Higgs doublet; for run = 350 GeV, it translates into a ~ 30% reduction in 
BR(h — > 77), already in some tension with data. Non-decoupling _F-term models, consisting of new 
chiral superfields with hypercharge zero, predict < 1. 

In Sec. 4 we give some natural SUSY predictions and compare them to data. In particular, we 
show that (i) significant enhancement of the h —¥ 77 rate, up to a factor ~ 4 times the SM value, is 
viable; (ii) the enhancement in h — > 77 cannot exceed the enhancement in h — > ZZ, WW by more than 
40%, in some tension with LHC data; (iii) the enhancement in VBF production cannot exceed the 
enhancement in GF production by more than 50%; (iv) the hint for ~ 100% enhancement in h — >• bb 
at the Tevatron is incompatible with natural SUSY. 

We conclude in Sec. 5. The appendix contains a discussion of a 125 GeV Higgs' implications for 
non-decoupling gauge extensions of the MSSM. 

2 Loop effects 

In this section we discuss the implications of light superpartners for Higgs production and decay 
channels at hadron colliders. We begin by summarizing the results. 

The effects we consider include loop contributions from the charged Higgs, higgsinos, gauginos, 
stops, sbottoms, and staus. Of these, the only quantitatively relevant effects (potentially larger than 
~ 5%) involve stops and charginos that affect the Higgs couplings to photons and gluons. The 
stop contribution to the Higgs-gluon coupling, r G , is directly related and opposite in sign to the stop 
contribution to the Higgs photon coupling, r* ; see Eq. (2.8). The chargino contribution to r 7 decouples 
for tan/3 > 5, as the chargino-Higgs coupling scales linearly with (l/tan/3). 

The stop and chargino effects are constrained by direct searches for these particles and, to some 
extent, by naturalness. Imposing collider limits together with demanding fine tuning no worse than 
1 : 10, we check the possible size of these effects by varying the relevant theory parameters. The 
results are plotted in Fig. 1, where the lines show the maximal and minimal values of r 7 obtained 
within natural SUSY. The possible stop contribution to rc is plotted in Fig. 2. 

The ranges of the different contributions to tq and r 1 in natural SUSY are summarized in Table. 1, 
limiting to tan/3 > 2. We also report upper limit estimates on the subdominant contribution of staus 
and sbottoms. In the following subsections we provide more details of the calculations. 
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Figure 1. The range in r 7 as a function of tan/3, compatible with fine-tuning no worse than 1:10 (details 
in Sees. 2.1.4 and 2.2). Left: stop contribution. Right: chargino contribution. In the left panel, solid curves 
correspond to total fine-tuning, defined in Eq. (2.13), while dashed curves correspond to tuning with respect 
to the Z boson mass alone, defined in Eq. (2.3). Note that the stop contribution to r 7 (left panel) is inversely 
related to re, see Eq. (2.8). 




Figure 2. Natural range in ra as a function of tan/3 (details in Sec. 2.1.4). Solid curves correspond to total 
fine-tuning, defined in Eq. (2.13), while dashed curves correspond to tuning with respect to the Z boson mass 
alone, defined in Eq. (2.3). Note that the stop contribution to re is inversely related to r 7 , see Eq. (2.8). 



2.1 Stops 

2.1.1 Higgs mass and fine-tuning 

The most widely studied radiative correction associated with stops is the increase in the Higgs mass. 
A well known, simple, analytic one-loop estimate that nevertheless includes the most relevant two-loop 
correction is given by [62] 
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Table 1. Ranges of re and r 7 in natural SUSY from different contributions, limiting to tan/3 > 2. The stop 
contributions to re and r 7 are anti-correlated, see Eq. (2.8). 



Here, X t — A t ~ ///tan/3, M 2 = mr rnj 3 , Qi = \/m t M s , Q2 = M s , and m t (Q) is the running top 
mass taken at scale Q, 

with M t being the top pole mass, taken here as 172.5 GeV. Eq. (2.1) slightly overestimates the Higgs 
mass, at a level of ~ 2 GeV, compared to the numerical two-loop package FeynHiggs [63]. Here in 
part of the analysis we will use Eq. (2.1) as a tight but robust upper bound to the Higgs mass in the 
MSSM, while in other parts we prefer the full FeynHiggs calculation, advising the reader accordingly. 

The correction to the Higgs quartic coupling, leading to Eq. (2.1), is accompanied by a shift to the 
Higgs self-energy that must be balanced by the bare mass in order not to destabilize the electroweak 
scale. Defining a fine-tuning measure, Az, with respect to the Z boson mass [64], we approximate the 
stop contribution to Az by 



(Az 1 ), 



25m% 



Sm 2 Hu \stop = ~^Vt (m 2 Q 3 + m 2 u3 + A 2 ) log (J^j) ■ ^ 

In what follows we assume conservatively a low SUSY breaking mediation scale A = 20 TeV. 

Eqs. (2.1) and (2.3) imply that a Higgs mass of 125 GeV would be highly fine-tuned in the 
MSSM. Unless we let go of natural SUSY, some physics beyond the MSSM must deform the scalar 
Higgs sector. Under this assumption, Eq. (2.1) does not directly constrain the stop sector. However, 
Eq. (2.3) continues to provide a reasonable guide to the fine-tuning associated with stops, with typically 
only mild modifications that we discuss in Sees. 3.3 and 3.4, in the context of concrete examples. 

2.1.2 Higgs couplings to gluons and photons 

The Higgs low energy theorem tells us that in the presence of heavy colored multiplets whose masses 
depend on the Higgs vacuum expectation value (VEV), e.g., the top-stop multiplet, the leading log 
Higgs-gluon-gluon coupling at the one- loop level is given by [65, 66] 



r - — 

'-hgg — ^2 



s _h f 2 y tF d\o g detm F (v) + ly ts dlogdet ms (v) \ 
TT V \ o log v 2 ^ o log V J 

where F and S denote respectively colored fermion and scalar with Dynkin index tx (= 1/2 for the 
fundamental representation) and mass matrices m(v). Applying Eq. (2.4) to the top-stop multiplet, 
we have [67] 
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where we neglect D-terms. We factored out an overall correction factor r t , defined via Eq. (1.1), that 
comes about by Higgs mixing. The total hGG vertex correction reads 

r G = r t r G . (2.6) 

Eq. (2.5) compares well with results from FeynHiggs; nevertheless, in numerical computations we 
include the £>-term contribution. Concerning the leading log approximation, this can be checked by 
comparing the full fermion and scalar loop function ratio evaluated at m t and mi, 
\Fq (m^/Am-) /J r i/2 ( TO h/4m 2 )], to the asymptotic value (1/4) of this ratio at — > 0. Varying 
between 150-1000 GeV, we find that the leading log approximation is good to about 6%. 

For reasonably light stops, Eq. (2.5) leads to a substantial effect, e.g. with m; = = 250 
GeV, tq = 1.24, implying 53% increase in GF production. As long as stop mixing is small, the hGG 
coupling is enhanced compared to the SM and consequently the GF rate is enhanced. As discussed 
in [66], large X t could in principle reduce the hGG coupling. However, naturalness, together with the 
direct bound mj > 100 GeV, limit this possibility: large X t adds to the weak-scale fine tuning both 
directly, through Eq. (2.3), and indirectly because it requires a larger diagonal soft mass to start with. 
In the next section we exhibit further constraints on such large X t that arise from rare B decays at 
large tan /3. 

There is an inverse correlation between the top/stop contributions to the Higgs effective coupling to 
photons and to gluons, the negative sign coming because of the dominant W diagram that contributes 
to /177 with opposite sign from the matter loops. To see this, let us denote the W and top loop 
contributions to the /177 amplitude by A n w and AJ , respectively, and the stop contribution by A~. 
Let us further define the hGG top and stop- induced amplitudes by Af and Af , and note that 

A] Af T 

^ = J?=rh-l, (2-7) 



to leading order in a s . This gives 
17 



Alir + Al + At 



(Al + A]f M 



1.28ry — 0.28r G , W, top, and stop contributions, (2-8) 



using A^ w « 8.33 and Al « -1.84 in the SM, valid for m h = 125 GeV. 

Eqs. (2.6) and (2.8) do not include loop contributions of additional particles, notably charginos 
and bottom and tau fermions and scalars. The bottom and tau fermion contributions remain below 
about five percent of the top even for rt>. T ~ 10. The chargino, sbottom and stau contributions can in 
principle become relevant in some corners of the MSSM parameter space, resulting with some loss of 
predictivity by disturbing the r 7 — tq correlation of Eq. (2.8). Below we examine these terms in more 
detail, concluding that in natural SUSY, the sbottom and stau contributions can be neglected while 
charginos may lead to marginally observable effects. 

2.1.3 Large stop mixing vs. fine-tuning in BR(B — » X s j) 

Light, mixed stops are constrained by rare B decays. The branching fraction for the rare decay 
B — > X s j has been measured experimentally to a precision of better than ten percent [68], 

BR(B ->■ V s 7) exp = (3.52 ± 0.25) x 10~ 4 . (2.9) 

The theoretical SM NNLO calculation has reached a similar accuracy [69] 2 , 

BR(B -> V s7 ) SM = (2.98 ± 0.26) x 10~ 4 . (2.10) 

2 Ref. [70] found the theoretical result BR(B -> X s -y) SM = (3.15 ± 0.25) X 10" 4 . 
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The theoretical NNLO SM prediction is fully determined by observable quantities, namely the masses 
of the top quark and W boson and gauge couplings. Therefore, the agreement (within ~ 1.5a) between 
Eqs. (2.9) and (2.10) allows us to define an observable quantity, O bsi , that we can compare against 
models of new physics, 



_ BR{B^X sl )^ 

~ BR(B^X sl )SM " 1 " 0-18 ± 0.13. (2.11) 



Eq. (2.11) means that new physics is now only allowed to contribute to B — > X s ^ at about 
thirty percent of the SM contribution. Because the SM contribution begins at one-loop, new physics 
models such as SUSY can easily produce larger contributions. Recalling the possibility of accidental 
cancellations, typical SUSY Higgs analyses in the literature either ignore B —> X s j or focus on 
parameter regions where cancellations occur. Here, given our interest in natural models, we will use 
Eq. (2.11) to estimate the level of fine-tuning involved in the latter approach [71]. 

Given a model input parameter P (e.g., A t ) that contributes to BR(B — > X s ^f), we assess the 
degree of fine-tuning A associated with it in a similar way to the fine-tuning measure commonly 
associated with the Z boson mass. The only slight modification we apply here is to account for the 
uncertainty in the experimental determination of 0f, S7 : 



A:; 1 

<->bs~f 



O bsi <91og0 6s7 




P dO bsi 


o<D bai d\ogP 




0.3 dP 



(2.12) 



where we chose to combine linearly the absolute values of the central value and of the uncertainty on 
the right hand side of Eq. (2.11), setting <jQ b — 0-3- A total fine tuning is defined as 



A-^/A^ + A^. (2.13) 



Consider now minimally flavor violating new physics contributions to the Wilson coefficients CV,8 
of the electromagnetic and chromomagnetic dipole operators 07,8 j 



°7 = T^ 771 ^ (*L<T^b R ) F^, 8 = ~^m b {s L a^T a b R ) G a ^ . (2.14) 
loir^ 1671"^ 



Taking 6*7,8 to be input at the top mass scale, the contribution to O bsi can be approximated by [72] 

O bsi = a 77 \C 7 \ 2 + a S8 \C 8 \ 2 +M{a 7 C 7 + a s Cs + a rs C 7 C* 8 }, (2.15) 

with a 7 = -2.41 + 0.2H, a 8 = -0.75 - 0.19i, a 77 = 1.59, a 88 = 0.26, a 78 = 0.82 - 0.30i. The 
MSSM with light mixed stops and light higgsinos makes a tan /3-enhanced contribution to the Wilson 
coefficients, of the form 



cv,8 ~ — —^7,8 r~l2>"T~l2 tan ^' 2 - 16 
mf \ \fi 2 \n 2 ' 



where F 7 $(x,y) are loop functions that take 0(1) values for x ~ y = 0(1). Using Eqs. (2.16), (2.15) 
and (2.12) we see that if we wish to avoid accidental cancellations to a level of one part per ten, 
then we must have (A t ^(tan/3/m?) < few. In our analysis we compute the stop contribution to 0f, S7 , 
and exclude a parameter space from being part of natural SUSY based on a tuning criterion, such as 
A to t < 10%, rather than on the precise deviation of O bsl from its central value. 



-7- 



2.1.4 Summary of stop effects 

In Figs. 1 (left) and 2 we study the range of values for r* and r G that can be naturally expected 
due to light stops. Here, we plot the extremal (upper and lower) values of r* and r G that can be 
achieved by varying A t , /x, mQ,mjj for different values of tan /3, while imposing A tot > 10%, or tuning 
no worse than one part per ten. We also, very conservatively, demand that > 100 GeV. The extra 
fine-tuning incurred by large A t , as seen by Eqs. (2.3) and (2.16), suggests that stop mixing should 
not be large, and more so for larger tan/3. As a result, the lower bound on r G (upper bound on r* ) is 
set by naturalness. 

We learn that r G < 0.85 does not arise within natural SUSY, but a rather significant increase of 
r l G ~ 1.4 is obtained with light, unmixed stops. The upper bound r G < 1.5 obtains when stops are un- 
mixed and as light as we let them, close to the top mass. Ignoring the possibility of additional SUSY 
contributions, this ultra-light region is in strong tension with electroweak precision tests (EWPTs). For 
example, setting X t = 0, toq = mjj — 100 GeV (mi w 200 GeV) and neglecting sbottom mixing we 
find (Ap/p) 3.5 x 10~ 3 , more than 8a deviation. However, it is not inconceivable that contributions 
from the Higgs and gaugino sectors could ameliorate this tension. While such cancellation would 
certainly be fine-tuned, the level of tuning does not quite make it to the 1:10 level except at this very 
low mass end. Instead of complicating the analysis by folding in EWPTs into our naturalness criterion, 
we simply report that if we impose that the stop-sbottom contribution to (Ap/p) remains within the 
4cr range, we obtain stronger limits on the Higgs-gluon and Higgs-photon vertices, 0.9 < r G < 1.3 and 
0.9 < r\ < 1.03. 

Finally, while our philosophy in this work is that additional physics beyond the MSSM must affect 
the Higgs sector to account for the Higgs mass, it is nevertheless interesting to exhibit the implications 
of naturalness for the stop contribution to mu- In Fig. 3 we plot the maximal value of rrih obtained 
by varying the SUSY parameters as above. We use the analytical estimate for rrih and so this plot 
provides an upper bound. The plot shows how tuning for B — > X s j makes a large stop contribution 
to rrih less plausible for large tan j3. Fig. 4 shows contours of the MSSM Higgs mass (computed now 

using FeynHiggs), the total fine tuning defined in Eq. (2.13) and (j^jj m Eq. (2.5), for two values 
of tan/3 = 10,30. At large tan/3, tuning for B — > X s f disfavors the high mixing region, causing the 
maximal value of rrih to drop substantially. 

2.2 Charginos 

Naturalness dictates that at least one chargino must be light, m^± < 200 GeV. Hence, the chargino 
contribution to r 7 may be expected to become relevant [73-75]. What limits the effect to be modest 
is the direct bound, that we take to be rn^± > 94 GeV [76]. Imposing this bound, we compute 
the chargino contribution to r 7 , varying the relevant parameters in the range —300 < p/GeV< 300, 
< M 2 /GeV< 1000, 1 < tan/3 < 40. The result is shown in the right panel of Fig. 1. We conclude 
that: 

• A sizable effect is possible only for low tan (3 < 3, mostly limited to a reduction in r 7 . The effect 
is largest for tan/3 = 1, where we find 0.7 < r* < 1.13. If we restrict to tan/3 > 2, we have 
0.8 < < 1.1. The sign of (r* — 1) depends upon the sign of (/X/W2). 

• Restricting to tan /3 > 3 (tan j3 > 5) diminishes the effect, as the chargino-Higgs coupling oc sin 2/3 
is reduced. Here we find |r* — 1| < 10%(6%), where saturating the upper limit requires two 
very light charginos with m~± <~ m~± ~ 100 GeV. 
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Figure 3. Upper limit on van in the MSSM, as function of tan/3. The solid line corresponds to scenarios with 
up to 10% total fine-tuning (defined in Eq. (2.13)), while the dashed line corresponds to Z boson mass tuning 
alone (Eq. (2.3)). 
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Figure 4. Contours of the Higgs mass, the total fine tuning and Vq in the (mQ 3 ,Xt) plane. We set n = 150 
GeV, rriQ 3 = m U3 and X t = A t — /i/tan/3. 



The upper bound of ~ +10% for the chargino enhancement to the Higgs-photon coupling can be 
understood as follows. Taking tan/3 = 1 one obtains, at leading log, 



-6.49 



8 



' 3 M 2t i- m j v 8 m w I 2m w \ 

-6.49 ~ 3 x 6.49 to 2 ± I m~± j ' { ' ' 



where we took M2H > m^y in order to obtain positive interference with the VF-dominated SM ampli- 
tude Asm ~ —6.49. Imposing the bound m~± > 94 GeV, we obtain — 1 < 10%, in good agreement 
with the full one-loop computation. 

The bottom line is that a chargino contribution to the /177 vertex, with a sign that is theoretically 
unconstrained, can disturb the correlation between r 7 and re that we have found, accounting for the 
stop contribution alone. As a result we will be forced to assess rc and r 7 independently when we 
come to predict Higgs observables. Nevertheless, it will still be useful to describe the modified /177 
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and hGG in terms of Tq and . This separation becomes practical for tan j3 > 3, when the chargino 
correction decouples, with ramifications to Higgs mixing effects. 

2.3 Charged Higgs 

A charged Higgs loop diagram contributes to r 7 , with [77] 

/ Kh+h- \ ( m H± 

1 \MSSM 
\ A hH+H~ 

where ^h+h- = (# 2 _ s' 2 )/ 4 ~ 0.07 is the MSSM coupling. The contribution is negligible unless the 
charged Higgs is very light, in strong tension with B — > X s ^ [70] 3 , or the coupling Xhh+h- receives 
very large corrections from an extended Higgs sector. We note that in i^-term models where singlet 
chiral superfields are added, to be discussed in more detail shortly, a numerically large correction to 
some Higgs quartic couplings is conceivable; however, the coupling X^h+h- remains unaffected. 

2.4 Staus and sbottoms 

The possibility that light scalar t could boost the h/yy coupling was entertained in [78, 79]. This 
possibility is, however, outside of the scope for natural SUSY. Naturalness limits this effect as follows. 
The stau mass eigenvalues are 



m 2 +m 2 + rnzoze 

< 2 = — o — ± 



(m 2 f -mj+ 2m 2 z s 2 v c 2 p) 
\\ ^ ^ —t- + rnlXl (2.19) 



with Xf = Af — fj, tan /3. Let us neglect D-tevuis in the following discussion, for clarity, including them 
numerically later to verify our conclusions. The leading log contribution to r 7 can then be estimated 
as 

r 7 - 1 « - d]o *(<0 = l + stau contribution , (2 . 20 ) 

6A S m alogv 3 Asm m i x m f 2 

where Asm ~ 6.49 is the SM (W and t loop) amplitude, evaluated at rrih = 125 GeV. Using the fact 
that m~ 2 > + 2m T \X?\, we obtain 

< r 7 - 1 < Ar T x 10^ ( J^L) ( 1Q ^)" , stau contribution. (2.21) 

Imposing conservatively > 100 GeV and demanding fine-tuning no worse than 5% < 300 
GeV), we find that in natural SUSY the effect is not larger than 3% for tan/3 < 50 and r T ~ 1. 

The sbottom contribution to r 7 is estimated to be even smaller because of the smaller electric 
charge, 

< r 7 - 1 < r b x 10- 4 ij^j do^y) ' . sbottom contribution. (2.22) 
The contribution to the hGG vertex is given by 

r G - 1 m - ( ' ' 2 | , sbottom contribution, (2.23) 



'The tension with B — > X s -f may be ameliorated if the charged state is taken from an additional inert Higgs multiplct. 
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that we can bound by 

> r G - 1 > -0.01(r 6 /r t ) ^ J^j,^ ( 10 ™g cV ) " ' sbottom contribution - ( 2 - 24 ) 

To make the stau or sbottom contributions to r 7 or tq larger than a few percent, large values of 
r T ,r b , arising from Higgs mixing, would be required. As we show in the next sections, (i) one expects 
fb ~ T r ; (ii) large r b would increase the total Higgs width, implying a suppression to the h — > 77 rate 
that would more than compensate for the presumed loop correction and that is not seen in the data; 
and (iii) r& > 3 does not arise in the MSSM or any extension we are aware of for lifting the Higgs 
mass. We conclude that the sbottom and stau contributions to r 7 and tq can be safely neglected in 
assessing the predictions of natural SUSY to a few percent accuracy. 



3 Higgs mixing 

Including quantum corrections, the natural Higgs mass prediction in the MSSM does not exceed 100- 
110 GeV (see Fig. 4). A Higgs at 125 GeV then implies that quartic couplings in the scalar potential 
receive corrections at the ~ 50% level from new physics beyond the MSSM. These corrections can 
modify the Higgs-fermion and Higgs-vector couplings from their MSSM values at a similar level, in 
a model dependent manner, through Higgs mixing. In this section we derive the Higgs coupling 
modifications due to Higgs mixing. We do not commit ourselves to the MSSM structure, making 
instead the more general assumption of a 2HDM framework. We then discuss the implications for 
specific model examples. 

3.1 Two Higgs Doublet Model 

We do not consider here the possibility of additional light fields mixing with H Uj( i, as could occur 
e.g. if the MSSM is augmented with a light singlet. In addition, we assume that non-renormalizable 
interactions in the effective 2HDM can be neglected. Then above the weak scale but below (at least 
most of) the superpartners, the supersymmetric Higgs sector is described by an approximately type-II 
2HDM [80, 81], 

- C = H\V 2 H 1 + h\v 2 H 2 + mllHil 2 + m 2 2 \H 2 \ 2 

+ y \Hi\ 4 + y |ff 2 | 4 + \ 3 \H x f\H 2 \ 2 + XilH^^l 2 
+ {^{H\H 2 ) 2 + {H\H 2 ) (m 2 12 + X^ 2 + \ 7 \H 2 \ 2 ) 

+ Y t H 2 ei R Q L3 + (Y b H\ - Y b A b H^ b R Q L3 + (y t H\ - Y T A T H^ f R L 3 + h.c.} , (3.1) 

where H\ 2 ~ (l,2) +1 / 2 and we identify H 2 = H u , Hi = io 2 H* d . In Eq. (3.1) we include "wrong" 
Higgs couplings for bottom and tau fermions with coefficients Ab jr , that come about from integrating 
out third generation squarks, higgsinos and gauginos at one-loop. We omit first and second generation 
fermions. 

Traditionally, the analysis of Higgs couplings is achieved by diagonalizing the Higgs mass matrix 
and expressing the couplings in terms of the rotation angle a, connecting the interaction basis of 
Eq. (3.1) to the mass basis, and the ratio [80] 

tan/3=-^4. (3.2) 
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Omitting for the moment A;, jT , we have r b = r T , and 



V9hbb _ sma v 9hti _ cos a vg hvv 

implying the inequalities 



»"& = = 3, r t = = -— r v = = sm (/3 - a) , (3.3) 

mb cos p nit sm p zto^ 



r\ < tan 2 /3 + 1, r\ < — + 1, < 1. (3.4) 
tan p 

We are free to choose two independent parameters to describe r b , r t and ry. We choose these 
parameters to be tan/3 and r b . With this choice we write 



valid for all tan/3. We assumed that r t > 0, taking the positive root. 

We now comment on the validity of neglecting A(,. T in Eq. (3.5). As it turns out, Ab jT , that enter 
the Higgs couplings tan/3-enhanced, only become quantitatively relevant at large tan/3. However, if 
(thj tan/3) 2 1, then deviations in rt, ry are suppressed compared to deviations in rf,. As we discuss 
below in more detail, this results in the fact that whenever the values of r t or ry are non-negligible 
phenomenologically, then Eq. (3.5) applies to high accuracy even when finite Ab. T are introduced. 
This conclusion is useful: it means that for arbitrary new physics deformations of the MSSM Higgs 
potential - just as long as the basic 2HDM structure is maintained - only two variables, and tan /3, 
are required to describe Higgs mixing effects on the lighter Higgs effective couplings. Note, finally, 
that similar diagrams to those that produce A& >T , also produce finite \j. As the latter couplings 
are vanishing in the MSSM, this can lead to non-negligible modification to Higgs couplings [77, 82]. 
However, these corrections are fully accounted for in Eqs. (3.3-3.5), by assuming renormalized couplings 
in the potential (3.1). 

For the purpose of understanding the phenomenology of specific models it is useful to express r b 
in terms of parameters in Eq. (3.1). As a simple but interesting scenario, consider tan/3 > 3, where 
we can use (l/tan/3) as an expansion parameter [77, 82]. We assume some hierarchy between the 
masses of the two doublets 4 , m 2 > m 2 ,, and neglect CP-violation and loop corrections from charged 
and pseudo-scalar Higgs states. Defining the quantities 

M? = ml + ^f, B = m? 2 + M, (3 . 6 ) 

where A35 = A3+A5 , a direct diagrammatic evaluation, treating the parameter (B/Mf) = (1/ tan/3) + 
0(1/ tan 2 (i) as perturbation, yields 



where we used m 2 H = (Mf ) + 0(1/ tan 2 0). The result for the coupling r T is similar to r^, replacing 
A fc -> A T . 

4 More precisely, we must require B 2 /(m? H — m^) 2 <C 1, where rriH.h are the mass eigenvalues. This condition can 

be put as tan /3 > (l - rn h/ m jj) 1 ■ 

5 Compared with the basis of [80], (B/M 2 ) ~ 1/ tan/3 and our A35 equals their A345. 
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The approximation Eq. (3.7) is useful because it allows us to understand generic predictions for a 
rather wide range of models [77, 82] . We will show in the next section that it is also fairly accurate. The 
couplings A5, A7 correspond to hard breaking of a U{\)pq symmetry under which (H u Hd) is charged. 
Many phenomenologically relevant models (e.g. the MSSM) break U(1)pq only softly and/or only at 
loop level, and thus A5J are suppressed while A3 can be 0(1). Hence in many cases, A35 controls the 
correction and its sign determines whether r& is increased or decreased from unity. Some generic 
classes of SUSY models predict fixed signs of A35 and thus the direction of the shift in TV 

In the next subsections we apply Eqs. (3.5) and (3.7) to SUSY models, including the MSSM and 
then extensions that can naturally accommodate rrih = 125 GeV. We show that even though 0(1) 
variations in r?, can easily occur in such models, nevertheless significant predictive power is maintained, 
particularly in the approximate PQ limit. 

3.2 MSSM analysis 

The tree-level MSSM quartic Higgs potential is given by G = SU(2)w x U(l)y D-terms, 
Vb = E f {h^TqHu + HlT§H d ) 2 = (K| 2 - |^| 2 ) 2 + • • • , 

G 

where the dots represent the charged Higgs potential, g and g' are the SU(2)w and the U(l)y gauge 
couplings. Mapping onto the 2HDM, we have at tree level 

A a = A 2 = -A35 = 9 + 4 9 « 0.14, A 5 = A 6 = A 7 = 0, A 6 = A T = 0. (3.8) 

The coupling A35 is negative, and so Eq. (3.7) tells us that the value of ri, is enhanced in the MSSM: 
r b w 1 + (m% + ml)/m 2 H w 1 + 0.25 (m# / 300 GeV) 2 . 

Finite Ab r and A5 ; 6,7 arise at loop order; using the results of [83, 84] we have, parametrically 6 




A 7 - -10~ 2 , A b - 3 x 10" 




A T ~ 10- 4 I ™ ) . (3.9) 



-3 

We can now get back to estimate the validity of some of our previous approximations. First, while 
Eq. (3.5) automatically accounts for finite As^j (either radiatively generated or otherwise), it does not 
encode finite A& and we should verify its aplicability. Note that A& enters Higgs effective couplings 
multiplied by tan/3. However, from Eq. (3.9) we expect that Eq. (3.5) should apply to reasonable 
accuracy even for tan /3 ~ 10, in which case, because of the relative tan 2 /3 suppression, deviations in 
r t , ry are already at most marginally relevant experimentally. In Fig. 5 we verify this point directly 
by comparing Eq. (3.5) with the full one-loop results from FeynHiggs. We find that for moderate 
tan/3 = 3 the agreement is excellent, while for tan/3 = 10 the agreement is better than 10%. This 
makes the overall error, incurred by using Eq. (3.5) for Higgs observables, no larger than order percent. 

Next, we consider the approximate PQ result, Eq. (3.7). As we argued in Sec. 2, naturalness 
favors At < 300 GeV. Thus unless tan/3 is very large, the contributions from A7 (and A^, A r ) are 
subdominant. In Fig. 6 we study the accuracy of Eq. (3.7) using the tree-level couplings (3.8) and 
comparing to the numerical results from FeynHiggs. We find that Eq. (3.7) captures the correct 
result to 5% accuracy for intermediate tan/3 = 10, and to better than 20% accuracy for large (small) 
tan j3 — 40 (3), where tan /3-enhanced loop corrections (terms 0(1/ tan 2 /3)) are in effect. 



6 A contribution to A7, coming from stau loops, may become comparable to the stop contribution if one allows very 
large (A-r/x/m?). 
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Figure 5. (rt — 1) as a function of tha for tan/3 = 3 (left) and tan/3 = 10 (right). In addition to the 
stop-higgsino loop with parameters At and fi as shown, a sbottom-gluino contribution is also included with 
m£>3 = 300 GeV, M g = 800 GeV. The solid curves are derived by using the value of r^, extracted from 
FeynHiggs, in our formula Eq. (3.5). The dashed curves are derived using FeynHiggs. Similar results are found 
for TV- 




Tan = 10, li = 150 GeV, m Q3 


= m„, - 300 GeV 


. A > = 


300 GeV (A„ = 13 %) 




200GeV(A„= 16%) 
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A,= 0(A,„, = 18 5 




Figure 6. rt — 1 as a function of tua for tan/3 = 40 (upper left), tan/3 = 10 (upper right), tan/3 = 3 (lower). 
The solid curves are tree-level leading results in the large tan /3 approximation obtained from Eq. (3.7) (setting 
A;, and A7 to be zero). The dashed curves are the full results computed with FeynHiggs. We decouple the 
right-handed sbottom in the computation. 



Last, we question the approximation ss r T , that is violated by 7^ A T . Again, the non- 
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holomorphic corrections are small unless tan/3 is large, so we can use Eq. (3.7). This gives 

n _ r bQ + (1 - r b0 )A;,tan/? 
r T no + (1 - r b0 )A T tan/3' 

with 7'bo = t t q denoting the result for = A T = 0. We conclude that « r T is a reasonable 
approximation as long as r^, r T = 0(1). The approximation can fail if the hbb and hrf couplings are 
strongly suppressed compared to the SM prediction; however, in that case we do not expect conclusive 
experimental information on these channels to become available any time soon. 

3.3 Non-decoupling D-term models 

In this section we study models with new gauge interactions, in which Z?-terms contribute the leading 
effect in raising the Higgs quartic couplings. Examples in the literature include [15, 85, 86]. We focus 
here on models in which the Higgs fields transform in a vector representation (H Ul H^) under the new 
gauge group, so that a fi term, fj,H u H,j, is allowed in the superpotential; we show that these models 
gcnerically predict an enhancement in the hbb and hrf couplings. Models that go beyond the vectorial 
charge assignment for the Higgs fields must gauge the U(1)pq symmetry and tend to consist a hybrid 
of Z?-term and F-term models [87, 88] that we discuss in the following section. 

Consider two-site Moose models with a product gauge group SU(N)a x SU(N)b and gauge 
couplings gA and gs (for example, SU{N) = SU(2) in [85]), depicted in Fig. 7. The product group 
is broken to the diagonal, identified as the SM electroweak gauge group G. This can be done, for 




Figure 7. The fundamental building block in vector-like D-term models. 

instance, by introducing bi-fundamental link fields (X, X), developing VEVs (X) = (X) through a 
superpotential W = A(detX + detX) + /TrXX. The broken generator vector multiplets acquire mass 
My- To arrange for a non-decoupling effect, large soft SUSY breaking mass M s can be introduced for 
the link fields. Then, the low energy effective potential of these models takes the same form as the 
MSSM D-term. potential with rescaled coefficients: 

- e f (- + «to*Y ^ a2(1+A, r (1+A,) (ki 2 - miT • 

(3.11) 

Further details of the derivation leading to Eq. (3.11) can be found in App. A. 
The effect on the Higgs quartic potential is simply: 

- Ai = -A 2 = A 35 = A-™ + £±±£pj . (3.12) 



- 15 - 



The key point is that the sign of A35 is always negative, as long as H u ,Hd are assigned charges in a 
vector-like representation. Therefore, by Eq. (3.7), these models generically predict an enhancement 
of rt,. In particular, for light stops and tan/3 > few, we can approximate A35W 2 rs — A|w 2 ~ ~ m \- In 
the absence of hard PQ-brcaking, A7 -C (A35/ tan/3), Eq. (3.7) is then simplified to 




(3.13) 



More specifically, accepting nih — 125 GeV, we can compute the size of the required A. In Fig. 8 
we set A = A' and plot contours of nth = 125 GeV for different values of A in the stop soft mass-mixing 
plane, using Eq. (2.1) modified by the Z?-term correction. Requiring light stops with fine-tuning no 
worse than 1:10, we find: 

A>0.5, (3.14) 

giving a deviation, (r& — 1), smaller by 10% than the value obtained from Eq. (3.13). As we comment 
in App. A, Eq. (3.14) implies that D-term models that address mh — 125 GeV should indeed be close 
to the non-decoupling regime. 



Contours of m h =1 25 GeV; tan|i=5 Contours of rrL=1 25 GeV; tan|i=20 




-600 -400 -200 200 400 600 -600 -400 -200 200 400 600 

X t [GeV] X ( [GeV] 



Figure 8. Contours of mi = 125 GeV in the (Xt,mQ 3 ) plane, for different values of the D-term correction 
parameter A. Left: tan/3 = 5; Right: tan/3 = 20. We have set A = A' for simplicity, but the results are not 
sensitive to this choice. 

Note that the D-term corrections lead to a slight numerical modification in our estimate of the Z 
mass fine-tuning, where, roughly, we should replace mz — > (1 + A)m^ in Eq. (2.3). Compared to the 
usual MSSM estimate, this allows for somewhat heavier stops, ~ 500 GeV, to still be consistent 
with tuning of order 1:10. This has a small effect on our results concerning the viable range for stop 
loop effects in /177 and hGG, allowing a slightly more pronounced decrease in r G (increase in r*). 

3.4 Non-decoupling F-term models 

Here we consider models that raise the Higgs mass through new interactions in the superpotential. 
The classic example includes a SM singlet, interacting with the Higgs doublets via 

SW = XSH u H d . (3.15) 
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(We will get the same results by making S a hypercharge-neutral SU(2) triplet.) If S is given a large 
soft SUSY breaking mass, m 2 ^> m? Hl M^, with tuh being the mass of the heavier MSSM Higgs doublet 
and M s a possible supersymmetric mass for S, then the effective potential below m s is modified with 
a non- decoupling correction, 

V = V MSSM + \X\ 2 \H u H d \ 2 . (3.16) 
This gives, in our notation of Eq. (3.1), 

A4 = A MSSM_ |A| 2 5 A3g = A MSSM + |A| 2, ^ 

By Eq. (3.7), these models tend to decrease r&. To estimate the size of the effect, note that the 
correction to the Higgs mass, still neglecting mixing, is 

t™t = rn%(J^ym 2 2l3. (3.18) 

In Fig. 9 we plot the value of A required for rrih = 125 GeV, by adding Eq. (3.18) to Eq. (2.1). Mixing 
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Figure 9. Contours of m/j = 120, 125, 130 GeV, in the (tan /3, A) plane. Smooth lines: tree level; dashed lines: 
including stop correction with Xt = 0, rrif 1 = m^ 2 = 380 GeV. 

with the heavy singlet S reduces the lightest Higgs mass due to level splitting [60]; as a result, for fixed 
m/j, the value of A in Fig. 9 serves only as a lower bound. We conclude that A35 is always positive in 
this model and much larger than its gauge-coupling value in the MSSM. This reduces the hbb coupling 
below its SM value. 

As pointed out in [60], the non- decoupling limit discussed above has limited applicability be- 
cause naturalness constrains m s <TeV. In contrast with D-term models, however, where electroweak 
precision tests constrain the SUSY scale, the singlet F-term example is phenomenologicallly viable 
also in the SUSY limit. In analogy with the D-term example, adding a supersymmetric mass term, 
6W D (M S /2)S 2 , the shift in A 35 is suppressed by factors of (m a /M s ); however, a supersymmetric cor- 
rection A7 = — (A 2 /x*/M s ) is generated. The A7 term modifies the Higgs mass by 5m 2 oc — (X7V 2 / tan/3). 
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With \M S \ = 700 GeV, = 200 GeV, tan/3 = 2 and |A| = 0.7, the Higgs mass could be raised to 125 
GeV [89, 90] with just a little help from ~300 GeV stops. From Eq. (3.7), the small negative A7 will 
also act to decrease r&. 

We should stress, however, that it is not difficult to construct F-term models that do not decrease 
rft. In particular, models that attempt to produce the /x term dynamically via a weak-scale singlet 
vacuum expectation value, tend to predict a light singlet, in which case the 2HDM analysis ceases to 
apply. An example is the Z3-NMSSM, 

W = XSH u H d + ^S\ (3.19) 

Parameteric scans of this model show that the shift in does not have a definite direction once S is 
allowed to be light [54, 91]. To understand how this can happen, consider Eq. (3.19) with the following 
hierarchy of masses: m 2 H « m% d > l^sl 2 > \ m H u \ 2 , where m 2 a < 0. Integrating out first H4 and then 
S, and expanding to leading order in (l/tan/3), we have 

[A[V v 2 { g*+g* [A[n 

6 \m 2 s \ + m 2 H -m 2 h \ 4 + 2\k\ 2 J ' [ ' 

Notice that there is no longer a correction scaling as (\X\ 2 v 2 /mjj); instead there are two opposite-sign 
contributions that could be parametrically comparable. 

Another example that increases the Higgs mass with no definite effect in the hbb coupling is found 
by adding an SU(2) triplet chiral superfield, A_ ~ (1, 3)_i, with super potential <5W = \A-H U H U . 
With a large SUSY breaking mass m s _, integrating out the field A_ inserts a hard SUSY-breaking 
correction A2 — > A2 + |A| 2 without modifying any of the other quartics in Eq. (3.1), lifting the Higgs 
mass with no further effect on the Higgs-fermion couplings. 

Our main conclusions from the discussion in Sees. 3.3 and 3.4 are these: (i) 0(1) modifications to 
r& )T are plausible through Higgs mixing in concrete extensions of the MSSM, that address the Higgs 
mass; (ii) measuring a deviation in r b will have strong implications for generic classes of models, in 
particular, non-decoupling D-term and _F-term models, that predict opposite sign effects. The value 
of (fb — 1) will provide very suggestive hints for the masses, quantum numbers and couplings of new 
particles beyond the MSSM. 



4 Natural SUSY predicts 

A 125 GeV SM-likc Higgs has several accessible decay and production modes. However, at this stage 
experimental uncertainties at the LHC and Tevatron are large. It is reasonable to estimate that 
even with the LHC 14 TeV run well under way (Ant = 30 ftV 1 ), individual couplings will only be 
measured to 20% — 40% accuracy [92]. These estimates are naive as far as systematics are concerned, 
but they should fall in the right ballpark. Given the experimental prospects, we limit the discussion 
to theoretical predictions of Higgs couplings that are valid to about 5%. This translates to ~ 10% 
accuracy in production and decay rates. Natural SUSY is perfectly capable of inducing much larger 
deviations, in which case our analysis will help to discriminate between different models. 

Following the discussion in Sees. 2 and 3, we classify the contributions to modified Higgs couplings 
into loop effects and mixing. We can now parametrize Higgs observables using four independent 
variables, two for each class of effects. These variables are 

tan/3, r b , (o.85 < r% < 1.5) , (o.7 < < l.lj . (4.1) 
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Using rfc and tan/3 in Eq. (3.5) we can compute r t and ry. Using r Gl , r t and ry we can compute 
the observable factors tq and r 7 . A set of four free parameters, with a limited range of values, is a 
rather predictive framework considering that experimental Higgs analyses will be sensitive to O(10) 
different production/decay channels. A few comments arc in order: 

• The fact that four variables suffice to describe Higgs production and decay is not special to 
our SUSY framework, but simply the result of assuming a 2HDM at the weak scale. This 
assumption gave us r t and ry in terms of and tan/3, 7 while to describe the Higgs photon and 
gluon couplings we could have chosen to use r 7 and tq directly. 



-± 



• Our choice of the r* and r G variables, is based on our ability to predict the viable numerical 
ranges for them in the particular framework of natural SUSY. It is of interest to spell out what 
part of the constraints on and r l G actually comes from naturalness vs. experimental limits: 

— The upper limit r G < 1.5 comes from imposing the direct constraint m 4 ~ >100 GeV and, 
when it is saturated, stops contribute significantly to (Ap/p). Requiring further that the 
stop-sbottom contribution to (Ap/p) does not exceed 4er would lower this bound to r G < 1.3. 
Similarly, both the upper and lower limits on the chargino contribution do not involve 
naturalness considerations, but merely the direct constraint m^± > 94 GeV. 

— The lower limit r G > 0.85 does arise from naturalness considerations. More specifically, 
it comes about by limiting the stop mixing to be modest. In some models, e.g. a singlet 
extension i^-term model with sizable A, the naturalness constraint may be relaxed somewhat 
and with it the lower bound on r G . 

— Finally, it was essentially naturalness (though assisted by direct constraints) that guided 
us to neglect the stau and sbottom loop corrections to r 7 and tq- 

Let us apply our analysis to a number of experimental channels, defining the signal strength 
p x = ° x BR(X) /SM. First, a 125 GeV SM Higgs has a partial width of w 64.4% to bb and rf, 24.3% 
to WW and ZZ, 8.5% to gluon pairs and 2.7% to cc. For our purpose it suffices to approximate the 
total width modification by 

r/r SM = Ptot « 0.64r 2 + 0.24r;> + 0.09r| + 0.03r t 2 . (4.2) 

Note that ptot depends mostly on Higgs mixing through and ry(r&, tan /3). Consider now the 
following six processes, with GF, VBF and AP standing for gluon fusion, vector boson fusion and 
associated production, respectively: 



2 2 



^T,GF = (4.3) 

2 2 

V lr ,VBF = -^f , (4.4) 

fJ-zz,ww;GF = r -f^~ , (4.5) 

^b.rr;AP = . (4.6) 

Some relevant questions, motivated in part by the current experimental situation (summarized in 
Sec. 4.1), are the following: 



7 Alternatively, of course, we could have used the angle variables j3, a. 
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1. What is the maximal enhancement for h — > 77 in GF production? The bound is obtained 
with large positive stop corrections, suppressed hbb coupling rj«l and SM-like r*t ~ Ty « 1. 
We find 

Vjt,gf < 4.7 to 6.2, (4.7) 

where we allow stops to provide r G =1.3 to 1.5, respectively. fij-y-GF ~ 2 is easy to achieve in 
beyond-MSSM SUSY by suppressed hbb coupling, together with a little help from natural stops. 
Note that our upper limits on the chargino loop contribution implies that if the diphoton rate 
is increased by more than ~ 20%, then so should the ZZ, WW rates. (See bullet (3) below). 

2. What is the maximal enhancement for h —> bb in AP? This question is partially motivated 
by the hint for factor ~ 2 enhancement in h — » bb associated production, reported recently by 
the Tevatron experiments. A constraint comes simply from 2HDM trigonometry, 

2 2 

Vbb;AP = ^ < 1-5. (4.8) 

Mtot 

We obtain the upper bound of 1.5 by examining 0.1 < tan f3 < 50, < n, < 10, and varying r G in 
the range of Eq. (4.1). The maximal value obtains for tan j3 > 30 and 77 > 3.5. Experimentally 
establishing a nontrivial lower bound on Hbb-AP would have profound implications for natural 
SUSY. For instance, establishing fibb;AP > 1-4 would provide the lower bounds tan/3 > 8, r-f, > 2, 
ruling out F-term models like ASUSY, while making a strong case for D-term models that 
enhance hbb. Interpreted within the latter models, by Eq. (3.13) the measurement would imply 
a challenging (but still consistent with direct searches) upper bound on the heavier Higgs doublet, 
niH < 230 GeV. This discussion applies also to /j, tt -ap- 

3. What is the maximal ratio of h — > 77 vs. ZZ, WW in the GF production channel? 

Both LHC experiments report /i 77; GF > 3 x ^ww-,gf, and a low WW rate is also found at the 
Tevatron. The ratio and the answer are given by 



M77;GF 



,.,2 



^WW,ZZ;GF r v 



1.28 - 0.28(r G r t /r v ) + (5r* fry) < 1.4, (4.9) 



where 6r* = — 1. The numerical value of 1.4 answers a slightly modified question, namely: 
"what is the maximal value of Eq. (4.9), assuming that h — > 77 is not suppressed, 
i.e. assuming /i 77: GF ^ 0.75?" This is more relevant, because to truly maximize Eq. (4.9) one 
would need to take ry small, which would diminish the actual observed 77 rate. Imposing 
fJ>-yj;GF > 0-75 and varying the variables in Eq. (4.1), we find the numerical bound in Eq. (4.9). 

4. What is the maximal ratio of h — > 77 in the VBF vs. the GF channels? Again, we 
further impose /i 77; gf > 0.75, this time in order to avoid irrelevant solutions with rj < 1 and 
vanishing GF production. In CMS, the reported VBF to GF ratio is ~ 2. The upper bound we 
obtain is: 

^■• VBF = T \< 1.5. (4.10) 

^77;GF r G 

The sample questions above give a sense to the level of predictive power in natural SUSY. However, 
the usefulness of reducing the number of free parameters to just a few comes mostly in interpreting 
a larger set of independent measurements. In what follows we demonstrate this point by fitting the 
parameters in Eq. (4.1) to current data and interpreting the results. 
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4.1 Interpreting current data 

Here we inspect and interpret current results from the LHC and the Tevatron experiments. We consider 
the following channels: 



1. 


CMS 


[26] 


and ATLAS [29] 77, including the CMS dijet tagged category. 


2. 


CMS 


[27] 


and ATLAS [30] ZZ. 


3. 


CMS 


[93] 


and ATLAS [94] WW. 


4. 


CMS 


[95] 


and ATLAS [96] rr. 


5. 


CMS 


[97] 


and ATLAS [98] bb. 


6. 


CDF 


[99] 


and DO [100] combined WW and bb. 



We included only the associated production mode for the rr and bb channels and the 77 dijet tagged 
category of CMS was interpreted as containing a 30% of events coming from GF and the remaining 
70% from VBF production. This is justified by the CMS estimate of the GF contamination, that can 
be found in [26]. The rates we use in this study were reported by the experiments as best fits to the 
signal strength and can be found in the Higgs combination papers and conference notes [31, 101, 102]. 

Our approach to interpreting the current data is as follows. First, we adopt a reference channel 
and assume some particular value for it, consistent with the experimental result. Then, assuming this 
reference value, we vary the parameters in Eq. (4.1) and obtain the viable ranges for all of the other 
experimental channels. We assume nih = 125 GeV, but our results have little sensitivity to varying 
m h by ±2 GeV. 

As the first reference channel, we select the LHC measurements of GF production h — »• 77. CMS 
and ATLAS report consistent results for this channel, with a relatively small experimental error, 
(Mt7;Gf) cxp ~ 1-5 ± 0.5. As reference value we choose: (fi^ Ti Qp ) rc = 1.5. The results are shown in 
Fig. 10. Within natural SUSY, our reference value for ^ 77; g_f is in about ler tension with each of the 
reduced WW and ZZ at CMS, the reduced WW at ATLAS, and the reduced (enhanced) WW (bb) 
at the Tevatron. Natural SUSY predicts that the ratio between the h — > 77 and h — > WW, ZZ rates 
should get closer to unity with future data. If this does not happen, then we will have an indication 
that some of our basic assumptions were wrong. A simple possibility would be e.g. that additional 
charged particles, beyond the MSSM matter content, contribute to the /177 vertex. 

Second, we select as reference the intriguing Tevatron hint for enhanced h — > bb, assuming 
(Hbb;ApY c — 1.5. The result is depicted in Fig. 11. In this scenario, natural SUSY predicts a 
strong suppression in the VV, 77 final state channels that cannot be compensated by MSSM particle 
loops, in tension with LHC data. 

5 Conclusions 

We present a detailed study of Higgs couplings in natural SUSY. Our approach is different than 
most existing analyses in that, while we do keep the discussion centered around the essential recent 
experimental inputs such as the favored value of the Higgs mass, we do not aim to fit the current 
data [92, 103-110]. Instead, our goal is to formulate predictions for Higgs observables in the context 
of a non-minimal Higgs sector that addresses the little hierarchy problem while staying in accord with 
m/j ~ 125 GeV. Our framework includes quantum corrections from MSSM particles and allows for an 
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Figure 10. Natural SUSY predictions, assuming GF production h — > 77 as indicated by the black crosses, 
are marked by thick bands. 
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Figure 11. Natural SUSY predictions, assuming associated production h — > bb as indicated by the black 
cross, are marked by thick bands. 



arbitrary 2HDM potential, without restricting to the MSSM structure that is likely distorted by new 
interactions. 

We separated the analysis into two types of effects: (i) loop corrections, affecting the couplings /177 
and hGG, and (ii) Higgs mixing, modifying the Higgs-fermion and massive vector boson couplings. 
In Sec. 2 we considered loop effects. Demanding fine-tuning no worse than about 1:10, we argued 
that only stops and charginos can induce coupling modifications larger than ~ 5%. Light, unmixed 
stops enhance hGG, the most relevant constraints being electroweak precision tests and direct stop 
searches. Demanding that (Ap/p) remains within 4a from the experimental value constrains the stop 
correction to r G < 1.3. Highly mixed stops can reduce the hGG coupling; here, bounds come from 
rare B decays and naturalness. Instead of imposing the bound due to BR(B — > X s -f) directly, we 
incorporate BR(B — > X s j) into our assessment of the model fine-tuning. This measure disfavors 
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light, mixed stops, with (X t/J, taxi P/rrv-) > 1. We find that r l G > 0.9 is preferred for naturalness. Light 
charginos can vary the hryy vertex in the range 0.7 < < 1.1, when both states lie close to the 
direct limit of ^100 GeV. The effect scales as oc (1/ tan/3) and goes below ±5% for tan/3 > 5. 

In Sec. 3 we considered Higgs mixing. We showed that once one assumes that the weak-scale Higgs 
sector is an approximately type-II 2HDM, with natural flavor conservation broken only by loop effects, 
then, to a good accuracy, the analysis of Higgs observables depends on only four free parameters, two 
describing Higgs mixing and two describing loop effects. Because of the large branching fraction of a 
125 GeV SM-like Higgs to bb, the hbb coupling r& is of key phenomenological importance. We estimated 
the size, and interpreted the implications of possible variations in rj,, expected within concrete SUSY 
models. Non-decoupling Z?-term models generically predict > 1, provided that the MSSM Higgs 
doublets are taken to transform in a vector representation of the new gauge group. The precise result 
depends on the mass of the heavier Higgs doublet and can be estimated as ~ 1 + 2(mf l /mH) 2 . For 
rriB = 350 GeV, consistent with all other experimental constraints, this gives r& ~ 1.25 that translates 
into a ~ 30% reduction in BR(h — > 77), already in some tension with current experimental results. 
Non-decoupling F-terxa models, consisting of new chiral superfields with hypercharge zero and ^TeV 
SUSY breaking mass, predict rj, < 1. 

In Sec. 4 we listed natural SUSY predictions to Higgs observables and interpreted the current 
data. In particular, we showed that: 

• Significant enhancement of the h — > 77 rate, up to a factor ~ 4 times the SM result, is viable; 

• The enhancement in h — > 77 cannot exceed the enhancement in h — > ZZ, WW by more than 
40%, unless the GF h — > 77 rate is itself reduced by more than 25% compared to the SM; 

• h — » bb in the associated production channel cannot be enhanced over the SM prediction by 
more than 50%. 

Our analysis captures the phenomenologically relevant Higgs couplings in natural SUSY models, using 
four free parameters with restricted numerical ranges. Considering that Higgs studies at the LHC and 
Tevatron yield O(10) experimentally independent observables, natural SUSY provides a predictive, 
falsifiable set-up, with promising opportunities to discriminate between model sub-classes in the near 
future. 

Note added: while this paper was being prepared for submission, a related analysis appeared [110] 
that partially overlaps with our discussion in Sees. 3.3 and 3.4. 
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A More details of D-term models 

Discussions of D-term models can be found e.g. in [86, 89, 111]. Here, for completeness, we repeat 
the derivation of the low-energy effective potential discussed in Sec. 3.3. Then we comment on the 
implications of Eq. (3.14) for this class of models. 
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The link fields (E, E) transform as bi-fvmdamentals under the product gauge group with a canonical 
Kahler potential 



(A.l) 



where qa,b are the gauge couplings of gauge group A and B and Va,b are the corresponding vector 
multiplets. After E, E develop a VEV, /, the product gauge group is broken down to the diagonal 
group, which is identified as the SM electroweak gauge group G with a massless vector multiplet 



V G 



gAV B + 9bVa 



V9 A + 9% 

The SM gauge coupling is given by g^ 2 = gj 2 + g^ 2 . The orthogonal combination, 

y _ -gAVA+ 9bVb 
y/9% + 9b 

is massive, with a mass term in the Kahler potential 

Ky = M 2 V 2 +■■■ , 



(A.2) 



(A.3) 



(A.4) 



where My = (g 2 A + g%)f 2 ■ Both MSSM Higgs fields are only charged under SU(N)a with an initially 
canonical Kahler potential 



K H =Y. Hje 9AVA Hi = £ (Hje^o Hi _9SM9A H f VHHi 



(A.5) 



where we expand to the leading order in Vj? . After integrating out Vh through its equation of motion, 
we have 



2 2 

9g9a 
a 9%M V 



(A.6) 



whose contributions to the Higgs quartic couplings are of order [i 2 /My. Electroweak precision tests 
require that My > 3 TeV [112]. Given the naturalness condition /i < 300 GeV, these supersymmetric 
corrections are thus too small to explain a Higgs mass at 125 GeV. 

Introducing SUSY breaking effects via a universal soft mass, M s , to the link fields, modifies the 
Kahler potential of the heavy vector multiplet, 



K v = {M v + 4 M 2 )V% 
Again after integrating out Vh, we have 



2 2 

9g9a 



1 



4 M 2 



K eff _ \p H^e 9aVG H - V 

h G A \K My+M* 

which gives rise to the low-energy potential in Eq. (3.11). 

Consider two limits of the modified Higgs couplings in Eq. (3.11): 



E H t T G H i 



(A.l) 



(A.8) 
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• My 3> M s , corresponding to a nearly supersymmetric heavy threshold, with 



MA g\ M 2 « g\ 



• M s ^> My, corresponding to hard SUSY breaking in the low energy potential, with 

A(A') = 4- 
9h 

Given that g^ 2 = g^ 2 +g^ 2 , we cannot take <?s very small at the weak scale. Assuming gA/gB = ^(1), 
we learn that JD-term models need a SUSY breaking mass M s at least comparable to the supersym- 
metric scale My. These models are thus unlikely to be well described by the nearly supersymmetric 
effective field theory of [89]. 

Finally we comment on the fine-tuning in these models. The new non-supersymmetric contribution 
to the Higgs quartic couplings leads to a quadratic divergence in the Higgs mass squared, cut off by 
M a . If we require this additional fine-tuning to be no worse than 10%, 

f|^4 < 10, (A.9) 
167T m A h 

we obtain a conservative upper bound, M s < 10 TeV. Taking into account the electroweak constraint, 
My > 3 TeV, and the requirement My < M s coming from Eq. (3.14), we find that M s is limited to 
the range 3-10 TeV. 
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